Search for K^K^ in ^p" Decays 
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KgK^l from tp" decays is searched for using the i/i" data collected by BESII at BEPC, the upper 
limit of the branching fraction is determined to be B{ip" K%Kl) < 2.1 x IQ-* at 90% C. L. The 
measurement is compared with the prediction of the 5*- and D-wave mixing model of the charmonia, 
based on the measurements of the branching fractions oi J/ip ^ KgKl and ip' K^K^. 
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I. INTRODUCTION dominated by the annihilation of cc into three gluons or 

a virtual photon, with widths proportional to the square 
of the wave function at the origin [3. This yields the 

From the perturbative QCD (pQCD), it is expected 
that both J/-0 and ip' decaying into light hadrons are 
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pQCD "12% rule", that is 

Qh = = -5 ~ l^/o- 

FoUowing the first observation of its violation in pi: 
and K*+K'^ + c.c. modes by Mark II 0, BBS has mea- 
sured many two-body modes of V'' decays, among which 
some obey the 12% rule while others violate it 0- There 
have been many theoretical efforts trying to solve the 
puzzle 0j, however, none explains all the existing exper- 
imental data satisfactorily and naturally. 

A most recent explanation of the "pTr puzzle" using the 
S- and D-wave charmonia mixing was proposed by Ros- 
ner In this scheme, the mixing of ■(/'(2'^S'i) state and 
il){l^Di) is in such a way which leads to almost complete 
cancellation of the decay amplitude of if)' pir, and en- 
hanced decay rate of ip" . A study on the measurement of 
ip" pn in e^e~ experiments shows that with the decay 
rate predicted by the S- and D-wave mixing, the interfer- 
ence between the three-gluon decay amplitude of the ip" 
and the continuum one-photon amplitude is destructive 
so the observed cross section is very small (jj, which is 
in agreement with the upper limit of the pn cross section 
at the ip" peak by Mark III 7]. Although this need to 
be further tested by high luminosity experiment operat- 
ing at the tp" mass energy, such as CLEOc H , it already 
implied that B{ip" pn) is most probably at the order 
of 10^'*, in agreement with the prediction of the S- and 
D-wsLve mixing scheme. 

If the S- and D-wave mixing is the key for solving 
the pTT puzzle, it applies to other decay modes as well, 
such as pseudoscalar pseudoscalar (PP) mode like KgK^. 
Recently, BES collaboration reported the branching frac- 
tions of KgKl final state in J/4> and ip' decays 9, lOj: 

BiJ/ip K^Kl) = (1.82 ± 0.04 ± 0.13) x lO"'', 
B{i}' K^Kl) = (5.24 ± 0.47 ± 0.48) x 10"^ 

These results yield Qk°k°^ — (28.8 ± 3.7)%, which is 
enhanced relative to the 12% rule by more than 4tT. 

By assuming that the pQCD 12% rule holds for K^Kl 
mode between J/^ and V'(2'^5'i), in the S- and D-wave 
charmonia mixing scheme, Ref. Ill| predicts the decay 
rate of ip" KgR^ in a range, that is 

0.12 ± 0.07 < 10^ X B{^" K%Kl) < 3.8 ± 1.1. 

Here the upper bound corresponds to (j) — 0° and the 
lower bound to (p — 180°, where (p is the relative phase 
between {K°Kl\l^Di) and {K'^Kl\2^Si). The uncer- 
tainties are due to the mixing ang le e between ip{2^Si) 
and ip{l^Di) states, and the measurements of B{ip' 
K%Kl) and BiJ/^p^K^sKl). 

In this paper, we report a search for ip" — > KgKl at 
BESII. 



II. THE EXPERIMENT 

The data used for the analysis are taken with the BE- 
SII detector at the BEPC storage ring in the vicinity 
of the Ip" peak (±4 MeV around ip" nominal mass). 
The data sample corresponds to a total of 17.7(1 ± 
5%) pb~^ luminosity as determined from large angle 
Bhabha events [T3 |. 

The BES is a conventional solenoidal magnet detector 
that is described in detail in Ref. .13il, BESII is the up- 
graded version of the BES detector [lj|- A 12-layer ver- 
tex chamber (VC) surrounding the beam pipe provides 
trigger information. A forty-layer main drift chamber 
(MDC), located radially outside the VC, provides tra- 
jectory and energy loss [dE/dx) information for charged 
tracks over 85% of the total solid angle. The momentum 
resolution is (Jp/p = 0.017-\/l -f {p in GeV/c), and the 
dE/dx resolution for hadron tracks is ~ 8%. An array of 
48 scintillation counters surrounding the MDC measures 
the time-of-flight (TOF) of charged tracks with a resolu- 
tion of ^ 200 ps for hadrons. Radially outside the TOF 
system is a 12 radiation length, lead-gas barrel shower 
counter (BSC). This measures the energies of electrons 
and photons over ^ 80% of the total solid angle with 
an energy resolution of ge/E = 22%/ Ve {E in GeV). 
Outside of the solenoidal coil, which provides a 0.4 Tesla 
magnetic field over the tracking volume, is an iron flux 
return that is instrumented with three double layers of 
counters that identify muons of momentum greater than 
0.5 GeV/c. 



III. MONTE CARLO 

Monte Carlo (MC) is used for mass resolution and 
detection efficiency determination, as well as the back- 
ground study. 

For the signal channel, ^" K^K^, the angular dis- 
tribution of or K'^ is generated as sin^ 9, where 9 is 
the polar angle in laboratory system. is allowed to 
decay according to its lifetime in the detector and only 
TT+TT^ is generated. For this study, 10000 events 
are generated. 

One of the main backgrounds is from e^e^ — > 77(7) 
events, with one photon converted into e^e~ pair in the 
detector material. This is studied with a MC sample 
which is 4 times as large as in the real data. 

Another background channel is K*^ {^92)K'^ -\- c.c. from 
ip" decays or e+e" direct production. 10000 events are 
generated for studying this background, which is about 
100 times more than in the real data. 

2M DD pairs are generated for studying the back- 
ground channels with c-quark production. This sample 
is about 14 times more than the real data sample. 
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The continuum channels from u, d, and s quark frag- 
mentation are generated with JETSET7.4 the MC 
sample is about 4 times as large as in the real data. 

The simulation of the detector response is a GeantS 
based package, where the interactions of the secondary 
particles with the detector material are simulated. Rea- 
sonable agreement between data and Monte Carlo simu- 
lation has been observed in various testing channels in- 
cluding e+e~ e+e~, e+e~ /^^M", J/4' ~^ PP and 
ij' TT+n-J/ij, J/^ ^ £+£- {£ = e, ^). 



IV. EVENT SELECTION 

The event selection criteria are all used in the analyses 
of the same final states at J/tp |^ and ip' energy [lol| . 
They are listed here for a easy reference. 

1. The number of charged tracks is required to be two 
with net charge zero. Each track should have good 
helix fit so that the error matrix of the track fit- 
ting is available for secondary vertex finding. The 
track is required to be within | cos6'| < 0.80, where 
9 is the polar angle of the track in MDC in the 
laboratory system. 

2. The two tracks are assumed to be 7r+ and tt~ , to 
find the intersect of the two tracks near the inter- 
action point, which will be taken as the secondary 
vertex. The tt'^tt" mass is required to be within 2a 
of the MC predicted mass resolution (8.9 MeV/c). 
The decay length in a;?;-plane, L^ij > 0.01 m is re- 
quired. 

3. The sum of the total energy of the photon can- 
didates < 1.0 GcV is used to remove the 
e~^e~ — > 77(7) backgrounds, with one photon con- 
verted into e"*"e~ pairs in the detector material. A 
neutral cluster is considered to be a photon candi- 
date when the angle between the nearest charged 
track and the cluster in the xy plane is greater than 
15°, the first hit is in the beginning 6 radiation 
lengths, and the angle between the cluster develop- 
ment direction in the BSC and the photon emission 
direction in xy plane is less than 37°. 

The above selection criteria are exactly the same as 
used for -0' KgK^ analysis To be unbiased with 
the expected small signal in the analysis, we try to fix 
the event selection criteria before looking at real data. 
By looking at the A'^ momentum distributions of the 
background MC samples generated, it is found th at t he 
background from the e+e" 77(7) and K*°{892)K" + 
c.c. is still very large. Then the selection criteria used for 
J/ip Ks^L selection are further applied to reduce 
the background level. 



4. The total BSC energy associated with the two 
charged tracks less than 1.0 GeV or the total XSE 
(the difference from the expected dE/dx for the 
electron hypothesis divided by the dE/dx resolu- 
tion) is less than —4. 

5. The opening angle between the two charged tracks 
larger than 20°. 

6. ElJ^ < 0.1 GeV, where ElJ^ is the sum of the ener- 
gies of the photon candidates outside a cone around 
the direction of (cos 6* < 0.95). 

Cuts 4 and 5 are used to reduce the gamma conversion 
background and cut 6 is used to reduce the K*^{%92)K'^+ 
c.c. background. After all the above cuts, there is no 
events left in the e+e^ 77(7) MC sample, and the 
remaining K*°{%%2)K° + c.c. events in the signal region 
is less than 1 after normalized to the A'*°(892)XO -I- c.c. 
cross section and the luminosity of the data sample. 

The signal region is defined as the momentum to 
be larger than 1.737 GeV/c, which is 2a {a = 42 MeV/c) 
lower than the MC predicted momentum for the sig- 
nal channel. 

After requiring all above criteria, the K'g momentum 
distribution of real data is shown in Fig. ^ as the black 
dots with error bars, the distribution of the events in the 
ATg mass sidebands is also shown (shaded histogram) , we 
can see that there is no clear difference between events 
in K'g mass region and those in K'g mass sidebands. In 
the same plot, we also give the MC predicted position of 
the signal events (blank histogram). It can be seen that 
there is no clear signal in data at the signal region. So 
we conclude that there is no K'gKl signal observed, and 
the upper limit of the tp" ^ ^s^l branching fraction 
will be determined based on the two candidates in the 
signal region. 



V. EFFICIENCIES AND SYSTEMATIC 
ERRORS 

The detection efficiency of the signal is estimated with 
10000 Monte Carlo simulated events, one gets emc = 
(41.44 ± 0.49)%, where the error is due to the statis- 
tics of the Monte Carlo sample. The trigger efficiency of 
KgKl events, which is lower because of the Kg decays, 
is measured to be (76.0 ± 1.8)% It is found that 

the reconstruction efficiency of the K'g in Monte Carlo 
is a bit higher than that in data |^], a correction factor 
of (96.3 ± 3.3)% should be applied to the Monte Carlo 
simulation. After taking into account all these factors, 
the global efficiency is e = 30.33%. 

The systematic error for the branching fraction mea- 
surement comes from the efficiencies of photon ID, 
secondary vertex finding, MDC tracking, trigger, the 



4 



2.5 



1.5 



0.5 







1.3 1.4 1 



5 1.6 1.7 1.8 1.9 
Kg momentum (GeV) 



2.1 



FIG. 1; The Kg momentum distribution. Data are shown by 
dots with error bars, the Kg mass sidebands background is 
shown by shaded histogram. The blank histogram is the MC 
simulated signal channel events, not normalized. 



VI. RESULTS AND DISCUSSION 



The upper limit of the branching fraction of 
K^Kl calculated with 



^obs 
''UL 



le 



N^nB{K% 7r+7r-)(l - a^v^) 



Where n^''£ is the upper limit of the observed number of 
events, which is 5.32 for 2 observed events at 90% C. L. 
assuming there is no background, e is the global efficiency 
for the signal channel. N^n is the number of -i/)" events, 
calculated with the total luminosity and the resonance 
parameters listed by PDG The systematic error of 
the measurement is considered by introducing 1 — u''^^ 
in the denominator of the formula for branching fraction 
calculation. 

Using numbers got above (listed in Table. ITT|l . one gets, 
at 90% C. L., 

B[il)" K%Kl) < 2.1 X 10"'*. 



branching fraction used, number of ip" events. Kg mass 
cut, angular distributions and so on. All these have 
been studied extensively in the analyses of J/-0 and 
i,' K^Kl H 113, and they are borrowed from the 
two analyses directly. 

Table. miists the systematic error from all sources. The 
uncertainties of the B{%Ij" e'^e^) and luminosity will 
affect the determination of the total number of the ip" 
events, the former comes from the PDG IT7| while the 
later is from the measurement used in Ref. Il2l. 



TABLE IL Numbers used in the calculation of upper limit of 
the branching fraction. 



quantity 


Value 




E 


5.32 
30.33% 






1.45 X 10^ 






0.6860 




sys 


17% 


B{i," ^ 


K%Kl) < 


2.1 X 10"* 



TABLE L Summary of the systematic errors. 



Source 


Systematic errors (%) 


MC statistics 


1.3 


Tptot 


2 




1.3 


2" vertex finding 


3.4 


MDC tracking 


4 


Trigger efficiency 


2.4 


B{xIj" ^ e+e-) 


15 


Luminosity 


5 


B{K% TT+TT-) 


0.4 


Total cr"^" 


17 



Add the errors from all the sources in quadrature, the 
total systematic error is 17%. 



Comparing with the corresponding theoretical calcu- 
lation of the branching fraction based on the S- and D- 
wave mixing model and the 12% rule, current upper limit 
is still well above the upper bound of the prediction 
of 3.8 X 10"^. To further pin down the upper limit of 
the branching fraction, thus to check the validity of the 
"Rosner's assumption" and the solution of the KgK^l en- 
hancement puzzle observed in ^' and J/tJj decays, a larger 
data sample is needed. The existing 55 pb~^ ip" data and 
the planned 3 Jh~^ tp" data samples from CLEOc Q will 
be obviously helpful for this study. 



VII. SUMMARY 

Flavor SU(3) breaking process K^K^ is searched for 
in tp" decays with BESll data sample at ip" energy, and 
the upper limit of the branching fraction is determined 
to be B{^" KlKl) < 2.1 x 10"^. The upper limit is 
still above the upper bound of the prediction if II] . 
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